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Abstract 

£ — . \ We explore the cosmological content of Salam-Sezgin six dimensional supergravity, and hnd a 

solution to the field equations in qualitative agreement with observation of distant supernovae, 
primordial nucleosynthesis abundances, and recent measurements of the cosmic microwave back- 
^ . ground. The carrier of the acceleration in the present de Sitter epoch is a quintessence field slowly 

rolling down its exponential potential. Intrinsic to this model is a second modulus which is au- 
tomatically stabilized and acts as a source of cold dark matter, with a mass proportional to an 
exponential function of the quintessence field (hence realizing VAMP models within a String con- 
text). However, any attempt to saturate the present cold dark matter component in this manner 
leads to unacceptable deviations from cosmological data - a numerical study reveals that this 
source can account for up to about 7% of the total cold dark matter budget. We also show that 
(1) the model will support a de Sitter energy in agreement with observation at the expense of a 
miniscule breaking of supersymmetry in the compact space; (2) variations in the fine structure 
constant are controlled by the stabilized modulus and are negligible; (3) "fifth" forces are carried 
by the stabilized modulus and are short range; (4) the long time behavior of the model in four 
dimensions is that of a Robertson- Walker universe with a constant expansion rate {w = —1/3). 
Finally, we present a String theory background by lifting our six dimensional cosmological solution 
to ten dimensions. 
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I. GENERAL IDEA 



The mechanism involved in generating a very small cosmological constant that satisfies 
't Hooft naturalness is one of the most pressing questions in contemporary physics. Re- 
cent observations of distant Type la supernovae [1] strongly indicate that the universe is 
expanding in an accelerating phase, with an effective de-Sitter (dS) constant H that nearly 
saturates the upper bound given by the present-day value of the Hubble constant, i.e., 
H < H ~ 1CT 33 eV. According to the Einstein field equations, H provides a measure of 
the scalar curvature of the space and is related to the vacuum energy density p vac through 
Friedmann's equation, 3 Mp X H 2 ~ p vac , where M P i ~ 2.4 x 10 18 GeV is the reduced Planck 
mass. However, the "natural" value of p vac coming from the zero-point energies of known 
elementary particles is found to be at least p vac ~ TeV 4 . Substitution of this value of p vac into 
Friedmann's equation yields H > 1CT 3 eV, grossly inconsistent with the set of supernova 
(SN) observations. The absence of a mechanism in agreement with 't Hooft naturalness 
criteria then centers on the following question: why is the vacuum energy needed by the 
Einstein field equations 120 orders of magnitude smaller than any "natural" cut-off scale in 
effective field theory of particle interactions, but not zero? 

Nowadays, the most popular framework which can address aspects of this question is 
the anthropic approach, in which the fundamental constants are not determined through 
fundamental reasons, but rather because such values are necessary for life (and hence intel- 
ligent observers to measure the constants) [2]. Of course, in order to implement this idea in 
a concrete physical theory, it is necessary to postulate a multiverse in which fundamental 
physical parameters can take different values. Recent investigations in String theory have 
applied a statistical approach to the enormous "landscape" of metastable vacua present in 
the theory [3]. A vast ensemble of metastable vacua with a small positive effective cosmo- 
logical constant that can accommodate the low energy effective field theory of the Standard 
Model (SM) have been found. Therefore, the idea of a string landscape has been used to 
proposed a concrete implementation of the anthropic principle. 

Nevertheless, the compactification of a String/M-theory background to a four dimen- 
sional solution undergoing accelerating expansion has proved to be exceedingly difficult. 
The obstruction to finding dS solutions in the low energy equations of String/M theory 
is well known and summarized in the no-go theorem of [4]. This theorem states that in 
a .D-dimensional theory of gravity, in which (a) the action is linear in the Ricci scalar 
curvature (6) the potential for the matter fields is non-positive and (c) the massless fields 
have positive defined kinetic terms, there are no (dynamical) compactifications of the form: 
ds 2 D = Vt 2 (y)(dx 2 d + g mn dy n dy m ) , if the d dimensional space has Minkowski 5*0(1, d— 1) or dS 
SO(l,d) isometries and its d dimensional gravitational constant is finite (i.e., the internal 
space has finite volume). The conclusions of the theorem can be circumvented if some of its 
hypotheses are not satisfied. Examples where the hypotheses can be relaxed exist: (?) one 
can find solutions in which not all of the internal dimensions are compact [5]; (ii) one may 
try to find a solution breaking Minkowski or de Sitter invariance [6]; (in) one may try to 
add negative tension matter (e.g., in the form of orientifold planes) [7]; (iv) one can even 
appeal to some intrincate String dynamics [8]. 

Salam-Sezgin six dimensional supergravity model [9] provides a specific example where 
the no-go theorem is not at work, because when their model is lifted to M theory the 
internal space is found to be non-compact [10]. The lower dimensional perspective of this, 
is that in six dimensions the potential can be positive. This model has perhaps attracted 
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the most attention because of the wide range of its phenomeno logical applications [11]. In 
this article we examine the cosmological implications of such a supergravity model during 
the epochs subsequent to primordial nucleosynthesis. We derive a solution of Einstein field 
equations which is in qualitative agreement with luminosity distance measurements of Type 
la supernovae [1], primordial nucleosynthesis abundances [12], data from the Sloan Digital 
Sky Survey (SDSS) [13], and the most recent measurements from the Wilkinson Microwave 
Anisotropy Probe (WMAP) satellite [14]. The observed acceleration of the universe is 
driven by the "dark energy" associated to a scalar field slowly rolling down its exponential 
potential (i.e., kinetic energy density < potential energy density = negative pressure) [15]. 
Very interestingly, the resulting cosmological model also predicts a cold dark matter (CDM) 
candidate. In analogy with the phenomenological proposal of [16], such a nonbaryonic matter 
interacts with the dark energy field and therefore the mass of the CDM particles evolves with 
the exponential dark energy potential. However, an attempt to saturate the present CDM 
component in this manner leads to gross deviations from present cosmological data. We 
will show that this type of CDM can account for up to about 7% of the total CDM budget. 
Generalizations of our scenario (using supergravities with more fields) might account for the 
rest. 



II. SALAM-SEZGIN COSMOLOGY 

We begin with the action of Salam-Sezgin six dimensional supergravity [9], setting to 
zero the fermionic terms in the background (of course fermionic excitations will arise from 
fluctuations), 

S=^-J d e x^\R - K 2 (d M a) 2 - K 2 e™F 2 MN - %e~™ - ^e 2 ™G 2 MNP ] . (1) 

Here, g 6 = det g M N, R is the Ricci scalar of g M N, F MN = d[ M A N] , G M np = d[ M B NP] + 
k,A[ M F NP }, and capital Latin indices run from to 5. A re-scaling of the constants: G 6 = 2k 2 , 
<f) = —kg and £ = Ag 2 leads to 

S = / d"x^[R - {d M <pf - ^ - ^e-4>F 2 MN - ^e^G 2 MNP ] . (2) 

The length dimensions of the fields are: [G 6 ] = L A , [£] = L 2 , [(f)] = [g 2 MN ] = 1, [A 2 M ] = L~ 4 , 
and [F MN ] = [G 2 MNP ] = L" 6 . 

Now, we consider a spontaneous compactification from six dimension to four dimension. 
To this end, we take the six dimensional manifold M to be a direct product of 4 Minkowski 
directions (hereafter denoted by Ni) and a compact orientable two dimensional manifold N2 
with constant curvature. Without loss of generality, we can set N2 to be a sphere S 2 , or a 
S 2 hyperbolic manifold with arbitrary genus. The metric on M locally takes the form 

ds 2 6 = ds 4 (t, x) 2 + e 2 ^W 2 , da 2 = { r l ffl + Sin ' f I l° T f (3) 
6 4V ' ' \ r 2 c (d$ 2 + smh 2 -&dip 2 ) for S 2 , w 

where (t, x) denotes a local coordinate system in A 1? r c is the compactification radius of N 2 . 
We assume that the scalar field <fi is only dependent on the point of N ± , i.e., = <f>(t, x). 
We further assume that the gauge field A M is excited on N 2 and is of the form 

r&cos* (s 2 ) 

^ - 1 b cosh d (E 2 ). 1 } 
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This is the monopole configuration detailed by Salam-Sezgin [9]. Since we set the Kalb- 
Ramond field B^p = and the term A^ M F NP ^ vanishes on N 2 , Gmnp = 0. The field 
strength becomes 

Fmn = Me-*f/rt . (5) 
Taking the variation of the gauge field Am in Eq. (2) we obtain the Maxwell equation 



d 



M 



V-4>p 



MN 



0. 



(6) 



It is easily seen that the field strengths in Eq. (5) satisfy Eq. (6). 
With this in mind, the Ricci scalar reduces to [17] 

R[M] = Rim] + e- 2 -fR[N 2 } - 4D/ - 6(<9 M /) 2 , 



(7) 



where R[M], R[Ni], and R[N 2 ] denote the Ricci scalars of the manifolds M, Aq, and N 2 ; 
respectively. (Greek indices run from to 3). The Ricci scalar of N 2 reads 



R[N 2 



+2/r 2 (S 2 ) 
-2/r c 2 (E 2 ). 



(8) 



To simplify the notation, from now on, Ri and R 2 indicate i2[Aq] and R[N 2 ], respectively. 
The determinant of the metric can be written as = e 2 ^ ' y/giyfg^, where g± = det g^ v and 
g a is the determinant of the metric of N 2 excluding the factor e 2 ^ . We define the gravitational 
constant in the four dimension as 
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(9) 



Hence, by using the field configuration given in Eq. (4) we can re-write the action in Eq. (2) 
as follows 



S = G~ I ^v^P^i + e~ 2f R 2 + 2(d,f) 2 - («9 M 0) 2 ] - J-e 2 ^ 



G ^ 2 e -2f-4> 



}. (10) 



Let us consider now a rescaling of the metric of Aq: g^ v = e 2 fg^ u and y/gl = ^fgl- Such a 
transformation brings the theory into the Einstein conformal frame where the action given 
in Eq. (10) takes the form 



1 

~G A 



d 4 xJg 4 R[g 4 ] - 4(<V) 2 - (d u 



-2f+4> 



G,b 2 



e -6f-* + e -*f R J. (11) 



The four dimensional Lagrangian is then 

R-4(d,f) 2 -(d„cf>) 2 -V(f,<f>) 



G 4 



(12) 



with 



V{f,<t>) = ^e- 2 ^ + 

(jTK 



G^_ e _ 6/ _ 



(13) 



where to simplify the notation we have defined: g = c/4 and R = R[gi\. 
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Let us now define a new orthogonal basis, X = (0 + 2f)/y/G A and Y = (<f> - 2f)/y/G 4 , 
so that the kinetic energy terms in the Lagrangian are both canonical, i.e., 



V9 



^r- l -{dxy- l -{dYf-v{x,Y) 



(14) 



where the potential V(X, Y) = V(f, 4>)/G± can be re-written (after some elementary algebra) 
as [18] 



V(X,Y) = 
The field equations are 



G,b 2 



IGlX 



G 6 



(15) 



R 



■\xv 



1 p _ G4 

2^" 2 



+ [ d ll Y(lY- f -^() 1l YO"Y) - tf /( „l (AM ) 



(16) 



□X = dxV, and DY = dyV. In order to allow for a dS era we assume that the metric takes 
the form 

ds 2 = -d1? + e 2h ®<B 2 , (17) 

and that X and Y depend only on the time coordinate, i.e., X = X(t) and Y = Y(t). Then 
the equations of motion for X and Y can be written as 



X + 3hX = -d x V 



and 



Y + 3hY = -dyV, 
whereas the only two independent components of Eq. (16) are 



h 2 
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and 



2h + 3h 2 ^^ 



■(X 2 + Y 2 ) +V(X,Y) 



--(X 2 + Y 2 ) + V(X,Y) 



(18) 
(19) 

(20) 
(21) 



The terms in the square brackets in Eq. (15) take the form of a quadratic function of 
e -V(hx^ This function has a global minimum at e~^^ x ° = R2 r^/(2 G$ b 2 ). Indeed, the 
necessary and sufficient condition for a minimum is that R2 > 0, so hereafter we only 
consider the spherical compactification, where e ~^^ x o — MpJ \Anb 2 ) . The condition for 
the potential to show a dS rather than an AdS or Minkowski phase is £b 2 > 1. Now, we 
expand Eq. (15) around the minimum, 



V(X,Y) 



s/GiY 

~g7~ 



!C + ^(X-X ) 2 + o((X-X f) 



where 



M x = 



t/k br c 



(22) 



(23) 
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and 

As shown by Salam-Sezgin [9] the requirements for preserving a fraction of supersymmetry 
(SUSY) in spherical compactifications to four dimension imply 6 2 £ = 1, corresponding to 
winding number n — ±1 for the monopole configuration. Consequently, a (Y-dependent) 
dS background can be obtained only through SUSY breaking. For now we will leave open 
the symmetry breaking mechanism and come back to this point after our phenomenological 
discussion. The Y-dependent physical mass of the Y-particles at any time is 

e VUlY/2 



M X (Y) = ^=^M X , (25) 

which makes this a varying mass particle (VAMP) model [16], although, in this case, the 
dependence on the quintessence field is fixed by the theory. The dS (vacuum) potential 
energy density is 

V Y = —^K. (26) 
In general, classical oscillations for the X particle will occur for 



M x >H=^^f^, (27) 

where p tot is the total energy density. (This condition is well known from axion cosmol- 
ogy [19]). A necessary condition for this to hold can be obtained by saturating p with Vy 
from Eq. (26) and making use of Eqs. (23) to (27), which leads to C,b 2 < 7. Of course, as 
we stray from the present into an era where the dS energy is not dominant, we must check 
at every step whether the inequality (27) holds. If the inequality is violated, the X-particle 
ceases to behave like CDM. 

In what follows, some combination of the parameters of the model will be determined by 
fitting present cosmological data. To this end we assume that SM fields are confined to Ni 
and we denote with p rad the radiation energy, with p x the matter energy associated with 
the Y-particles, and with p mat the remaining matter density. With this in mind, Eq. (19) 
can be re-written as 

Y + 3HY = (28) 
where V e s = Vy + p x and H is defined by the Friedmann equation 



H 2 = h 2 = 1 



3M 2 X 



-Y 2 + V eS + 

Prad fmat 



(29) 



(Note that the matter energy associated to the X particles is contained in V e g.) 

It is more convenient to consider the evolution in u = — ln(l + z), where z is the redshift 
parameter. As long as the oscillation condition is fulfilled, the VAMP CDM energy density 
is given in terms of the Y-particle number density n x [20] 

Px (Y,u) = M X (Y) n x {u) = C e^ Y ' 2 , (30) 
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where C is a constant to be determined by fitting to data. Along with Eq. (26), these define 
for us the effective (w-dependent) VAMP potential 



V cS (Y, u) = V Y + Px = A e^ Y + C e^ Y ' 2 e~ 3u , 



(31) 



where a A is just a constant given in terms of model parameters through Eqs. (22) and (24). 

Hereafter we adopt natural units, M P \ = 1. Denoting by a prime derivatives with respect 
to u, the equation of motion for Y becomes 



| 3y/ | d uP Y'/2 +3d Y V eS 



0, 



1 - Y' 2 /Q p 
where p = V e s + p ra d + p ma t- Quantities of importance are the dark energy density 

py = \h 2 y' 2 + v y , 

generally expressed in units of the critical density (Q = p/p c ) 

Py 



and the Hubble parameter 
The equation of state is 

Wy 



H 2 = 



3H 2 ' 
P 

3-Y' 2 /2 ' 



H 2 Y' 2 



Vs 



Y 



H 2 Y' 2 



Y 



(32) 

(33) 

(34) 
(35) 

(36) 



We pause to note that the exponential potential Vy ~ e Ay / Mpi , with A = y/2. Asymptotically, 



this represents the crossover situation with wy = — 1/3 [22], implying expansion at constant 
velocity. Nevertheless, we will find that there is a brief period encompassing the recent past 



< 



6) where there has been significant acceleration. 



Returning now to the quantitative analysis, we take p mat = Be~ Au and p rad = 
1CT 4 p mat e~ u f(u) [21] where B is a constant and f(u) parameterizes the M-dependent 
number of radiation degrees of freedom. In order to interpolate the various thresholds 
appearing prior to recombination (among others, QCD and electroweak) , we adopt a conve- 
nient phenomenological form f(u) = exp(— tt/15) [23]. We note at this point that solutions 
of Eq. (32) are independent by an overall normalization for the energy density. This is also 
true for the dimensionless quantities of interest fly and wy. 

With these forms for the energy densities, Eq. (32) can be integrated for various choices 
of A, B, and C, and initial conditions at u — —30. We take as initial condition Y (— 30) = 0. 
Because of the slow variation of Y over the range of u, changes in Y(— 30) are equivalent 
to altering the quantities A and C [24]. In accordance to equipartition arguments [24, 25] 
we take Y'{— 30) = 0.08. Because the Y evolution equation depends only on energy density 
ratios, and hence only on the ratios A : B : C of the previously introduced constants, we 
may, for the purposes of integration and without loss of generality, arbitrarily fix B and 
then scan the A and C parameter space for applicable solutions. In Fig. 1 we show a sample 
qualitative fit to the data. It has the property of allowing the maximum value of X-CDM 
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FIG. 1: The upper panel shows the evolution of Y as a function of u. Today corresponds to 
z = and for primordial nucleosynthesis z ~ 10 10 . We set the initial conditions Y{— 30) = and 
Y'{— 30) = 0.08; we take A : B : C = 11 : 0.3 : 0.1. The second panel shows the evolution of Vty 
(solid line), rj mat (dot-dashed line), and O ra d (dashed line) superposed over experimental best fits 
from SDSS and WMAP observations [13, 14]. The curves are not actual fits to the experimental 
data but are based on the particular choice of the Y evolution shown in the upper panel, which 
provides eyeball agreement with existing astrophysical observations. The lower panel shows the 
evolution of the equation of state wy superposed over the best fits to WMAP + SDSS data sets 
and WMAP + SNGold [14] . The solution of the field equations is consistent with the requirement 
from primordial nucleosynthesis, Qy < 0.045 (90%CL) [12], it also shows the established radiation 
and matter dominated epochs, and at the end shows an accelerated dS era. 



(about 7% of the total dark matter component) before the fits deviate unacceptably from 
data. 

It is worth pausing at this juncture to examine the consequences of this model for vari- 
ation in the fine structure constant and long range forces. Specifically, excitations of the 
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electromagnetic field on Aq will, through the presence of the dilaton factor in Eq. (2), seem- 
ingly induce variation in the electromagnetic fine structure constant a em = e 2 /47r, as well 
as a violation of the equivalence principle through a long range coupling of the dilaton to 
the electromagnetic component of the stress tensor. We now show that these effects are 
extremely negligible in the present model. First, it is easily seen using Eqs. (2) and (3) 
together with Eqs. (8)-(15), that the electromagnetic piece of the lagrangian as viewed from 
iVx is 

C em = —e-^fl , (37) 

where f^ v denotes a quantum fluctuation of the electromagnetic U(l) field. (Fluctuations of 
the U{1) background field are studied in the Appendix). At the equilibrium value X = X , 
the exponential factor is 

,— M 2 

~ Anb 2 ' W 

so that we can identify the electromagnetic coupling (1/e 2 ) ~ M^/b 2 . This shows that 
b ~ Mpi. We can then expand about the equilibrium point, and obtain an additional factor 
of (X — X )/Mpi. This will do two things [26]: (a) At the classical level, it will induce a 
variation of the electromagnetic coupling as X varies, with Aa cm /a cm ~ (X — X )/Mp\] 
(b) at the quantum level, exchange of X quanta will induce a new force through coupling to 
the electromagnetic component of matter. 

Item (b) is dangerous if the mass of the exchanged quanta are small, so that the force 
is long range. This is not the case in the present model: from Eq. (22) the X quanta have 
mass of 0(M x Mp\) ~ M P i/(r c 6), so that if r c is much less than 0(cm), the forces will play 
no role in the laboratory or cosmologically. 

As far as the variation of a cm is concerned, we find that px/Pm&t = (C / B)e Y /^ 2 ~ , so that 

^~3x 10- 120 e- 3w M^e y /^ 

= ±M* x (X-X ) 2 e Y ^M&. (39) 

This then gives, 



y/((X - X ) 2 ) = AX rms » 10- 60 e- 3 «/ 2 M P1 e y /( 2 ^/M x . (40) 

During the radiation era, Y ~ const ~ (see Fig. 1), so that during nucleosynthesis 
{u ~ —23) AX rms /Mpi ~ 10~ 45 /Mx, certainly no threat. It is interesting that such a small 
value can be understood as a result of inflation: from the equation of motion for the X field, 
it is simple to see that during a dS era with Hubble constant H, the amplitude AX rms is 
damped as e~ mt l 2 . For 50 e- foldings, this represents a damping of 10 32 . In order to make the 
numbers match (assuming a pre-inflation value AA rms /Mpi ~ 1) an additional damping of 
~ 10 13 is required from reheat temperature to primordial nucleosynthesis. With the e _3M ^ 2 
behavior, this implies a low reheat temperature, about 10 6 GeV. Otherwise, one may just 
assume an additional fine-tuning of the initial condition on X. 

As mentioned previously, the solutions of Eq. (32), as well as the quantities we are fitting 
to (fly and wy), depend only on the ratios of the energy densities. From the eyeball fit in 
Fig. 1 we have, up to a common constant, p rdinar y matter = Pmat oc 0.3 e -3 " and Vy oc 11 e v ^ 2Y . 
We can deduce from these relations that 

V Y( n0W ) = H e V2Y(now) _ 3 g e V2Y(now) _ 

p ma t(now) 0.3 
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Besides, we know that p mat (now) ~ 0.3p c (now) ~ 10 120 M^. Now, Eqs. (22) and (24) lead 
to 

W(now) = e^> ow ) (b 2 £ - 1) (42) 

so that from Eqs. (41) and (42) we obtain 

g^j <*-!)» 10-™. (43) 

It is apparent that this condition cannot be naturally accomplished by choosing large values 
of r c and/or b. There remains the possibility that SUSY breaking [27] or non-perturbative 
effects lead to an exponentially small deviation of 6 2 £ from unity, such that b 2 ^ = 1 + 
(9(10~ 119 ) [29]. Since a deviation of 6 2 £ from unity involves a breaking of supersymmetry, 
a small value for this dimensionless parameter, perhaps (1 TeV/M P1 ) 2 ~ 10~ 31 , can be 
expected on the basis of 't Hooft naturalness. It is the extent of the smallness, of course, 
which remains to be explained. 



III. THE STRING CONNECTION 



We now briefly comment on how the six dimensional solution derived above reads in 
String theory. To this end, we use the uplifting formulae developed by Cvetic, Gibbons and 
Pope [10]; we will denote with the subscript "cgp" the quantities of that paper and with 
"us" quantities in our paper. Let us more specifically look at Eq. (34) in Ref. [10], where 
the authors described the six dimensional Lagrangian they uplifted to Type I String theory. 
By simple inspection, we can see that the relation between their variables and fields with 

the ones we used in Eq. (2) is 0| cgp = -20| us , F 2 | cgp = VG^F 2 \ US , #3| cgP = \fG 6 /3G 3 \ us , and 
(8Gq)\ us . Our six dimensional background is determined by the (string frame) 



9 I cgp 
metric ds\ 



,2/ 



dt 2 



•2h 



dxi + r 2 da? 



the gauge field F, 



dip 



-frsin??, and the t- 



dependent functions h(t), /(f) = y/Ul(X - Y)/A, and <j>{t) = y/Gl(X + Y)/2. Identifying 
these expressions with those in Eqs. (47), (48) and (49) of Ref. [10] one obtains a full Type 
I or Type IIB configuration, consisting of a 3-form (denoted by F 3 ), 



8G 6 sinh p cosh p 
£ cosh 2 2p 
V2G 6 b 



dp A (da — \ TTTrb cos i3dip^j A (d(3 + J Tr^rb cos fidip) 



SG, 



^ cosh 2p 



sinfidO A dip A 



cosh 2 p \ da — 



8G, 



-b cos ddip 



sinh 2 p I d(3 + 



a dilaton (denoted by 0) 



8G ( 



-bcos$d(f 



2<t> 



2<j> 



cosh(2p) ' 

and a ten dimensional metric that in the string frame reads 



(44) 



(45) 



ds. 



str 



dsi + dz 2 + 



4G fi 



dp 2 + 



cosh 2 p 
cosh 2/5 



da — 



8G> 



-b cos fidip 
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') 



2" 



+ 



cosh 2p 



8G 6 



b cos fidtp 



(46) 



where p, z, a, and (3 denote the four extra coordinates. It is important to stress that though 
the uplifted procedure decribed above implies a non-compact internal manifold, the metric 
in Eq. (46) can be interpreted within the context of [7] (i.e., < p < L, with L > 1 an 
infrared cutoff where the spacetime smoothly closes up) to obtain a finite volume for the 
internal space and consequently a non-zero but tiny value for G 6 . 

IV. CONCLUSIONS 

We studied the six dimensional Salam-Sezgin model [9], where a solution of the form 
Minkowski 4 x S 2 is known to exist, with a £7(1) monopole serving as background in the two- 
sphere. This model circumvents the hypotheses of the no-go theorem [4] and then when lifted 
to String theory can show a dS phase. In this work we have allowed for time dependence 
of the six-dimensional moduli fields and metric (with a Robertson- Walker form). Time 
dependence in these fields vitiates invariance under the supersymmetry transformations. 
With these constructs, we have obtained the following results: 

(1) In terms of linear combinations of the S 2 moduli field and the six dimensional dilaton, 
the effective potential consists of (a) a pure exponential function of a quintessence field 
(this piece vanishes in the supersymmetric limit of the static theory) and (b) a part which 
is a source of cold dark matter, with a mass proportional to an exponential function of the 
quintessence field. This presence of a VAMP CDM candidate is inherent in the model. 

(2) If the monopole strength is precisely at the value prescribed by supersymmetry, the 
model is in gross disagreement with present cosmological data - there is no accelerative 
phase, and the contribution of energy from the quintessence field is purely kinetic. 
However, a miniscule deviation of 0(1CT 120 ) from this value permits a qualitative match 
with data. Contribution from the VAMP component to the matter energy density can be 
as large as about 7% without having negative impact on the fit. The emergence of a 
VAMP CDM candidate as a necessary companion of dark energy has been a surprising 
aspect of the present findings, and perhaps encouraging for future exploration of 
candidates which can assume a more prominent role in the CDM sector. 

(3) In our model, the exponential potential Vy ~ e Ay / M Pi ) with Y the quintessence field 
and A = \pl. The asymptotic behavior of the scale factor for exponential potentials 

e h(t) ^ go ^at for our case h ~ In t, leading to a conformally flat Robertson- Walker 

metric for large times. The deviation from constant velocity expansion into a brief 
accelerated phase in the neighborhood of our era makes the model phenomenologically 
viable. In the case that the supersymmetry condition (& 2 £ = 1) is imposed, and there is 
neither radiant energy nor dark matter except for the X contribution, we find for large 
times that the scale parameter e h ^ ~ so that even in this case the asymptotic metric 
is Robertson- Walker rather than Minkowski. Moreover, and rather intriguingly, the scale 
parameter is what one would find with radiation alone [28]. 

In sum, in spite of the shortcomings of the model (not a perfect fit, requirement of a tiny 
deviation from supersymmetric prescription for the monopole embedding), it has provided 
a stimulating new, and unifying, look at the dark energy and dark matter puzzles. 
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V. APPENDIX 

In this appendix we study the quantum fluctuations of the U(l) field associted to the 
background configuration. We start by considering fluctuations of the background field A° M 
in the 4 dimensional space, i.e, 

Am *■ A M + e a M , (47) 
where A M = if M ^ ip and a M = if M = ip. The fluctuations on A° M lead to 

Fmn F MN + e fMN ■ (48) 

Then, 

F MN F MN = g ML g NP [F MN F° LP + e F MN f LP + e 2 f MN f LP ] . (49) 

The second term vanishes and the first and third terms are nonzero because F MN ^ in 
the compact space and fuN ^ in the 4 dimensional space. If the Kalb-Ramond potential 
Bnm = 0, then the 3-form field strength can be written as 

Gmnp = kA[ M F N p] = — [Am Fjyp + Ap Fmn — An Fmp] ■ (50) 

Now we introduce notation of differential forms, in which the usual Maxwell field and field 
strenght read 

A 1 = A M dx M and F 2 = F MN dx M A dx N ; (51) 

respectively. (Note that dx M A dx N is antisymmetrized by definition.) With this in mind 
the 3-form reads 

G 3 = kA x AF 2 = kA m F np dx M A dx N A dx p . (52) 
Substituting Eqs. (47) and (48) into Eq. (52) we obtain 

G 3 = k [(A° M + ea M ) (F NP + ef NP ) dx M A dx N A dx p ] . (53) 
The background fields read 

A\ = b cos $ dip, F 2 ° = -b sin ■& d$ A dip , (54) 
and the fluctuations on the probe brane become 

d\ = a^dx^, j 2 = fdx^ A dx u , with / = d^a u — d^a u . (55) 

All in all, 
G 

—3 = A° V F^ d(p A d$ A dip + tA^f^ dip A dx^ A dx v + eF^a^ d§ A dip A dx» 

+ e 2 a fM f (v dx^ A dx ( A dx v . (56) 
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Using Eq. (54) and the antisymmetry of the wedge product, Eq. (56) can be re-written as 



K 



= e 



b cos fif^dip A dx^ A dx v - ba^ sin iW A dip A dx^ 1 + ea^f^dx* 1 A cfe c A dx u . (57) 



From the metric 



we can write the vielbeins 



ds 2 = e 2a dxl + e 2 \d$ 2 + sintf 2 V) 



(58) 



e a = e a dx a , = (Pd&, e v = e 13 sin ddip, 



dx a = e- a e a , d$ = e _/ V, dip 



sin-$ 



(59) 



where (3 = / + lnr c . (Lower latin indeces from the beginning of the alphabet indicate coordi- 
nates associted to the four dimensional Minkowski spacetime with metric r] ab .) Substituting 
into Eq. (57) we obtain 



G3 r, cos •§ 



K 



b -e- 2a -Pf ab e* A e a A e b - be' a ' 2 ^a a e^ A e v A e a + ee^ a a a f cb e a Ae c A e b , (60) 

sinv J 



where f ab = d a a b — d b a a . Because the three terms are orthogonal to each other straightfor- 
ward calculation leads to 



G\ = n 2 e 2 (b 2 cot 2 $e- Aa ~ 2 ?f 2 ab + b 2 e- 2a -^a 2 a ) + 0(e A ) . 
Then, the 5th term in Eq. (2) can be written as 

_L [ d 4 a .^ e 4a+2^ v ^ e -2« [ d $d V Sm$\U 2 t 2 b 2 COt 2 ^e'^) /, 

2Gfi J 6 J lv ' 



(61) 



Sg 3 - 



2 

ah 



+ (K 2 e 2 b 2 e- 2a -^)a 2 a 



(62) 



whereas the contribution from the 4th term in Eq. (2) can be computed from Eq. (49) 
yielding 



S F2 — 



J d 4 x^2ne 2 ^G 6 e 2 f ab 



2G 6 

-jd^ne 2 ^r 2 c e 2 f 2 ab 



Thus, 



Sg 3 + S F2 — — J 



d 4 x 



— f 2 + —a 2 

4 g2 J ab ' 2 



where the four dimensional effective coupling and the effective mass are of the form 



K 2 h 2 e -2<t> J Md(p gin ^ CQt 2 ^ 



OO 



and 



m 2 = 1™ We 2 "" 5 *- 2 * . 



(63) 



(64) 



(65) 



(66) 
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For the moment we let / dddtp sin d cot 2 d = N, where eventually we set N — > oo. Now 
to make quantum particle identification and coupling, we carry out the transformation 
a a — > gd a [30] . This implies that the second term in the right hand side of Eq. (64) vanishes, 
yielding 

fab = d a (ga b ) - d b (ga a ) = d a g a b - d b g a a + g d a a b - g d b a a = gf ab + a A dg (67) 
and consequently to leading order in N 

j 2 fl = j 2 tfft, + (a A dg) 2 + 2 g a b f ab d a g] . (68) 

If the coupling depends only on the time variable, 

^/aW a U (1)^ + 2^/* (69) 

where g = d t g and lower latin indices from the middle of the alphabet refer to the brane 
space-like dimensions. If we choose a time-like gauge in which a t = 0, then the term 
{g/g)dif u can be written as (l/2)(g/g)(d/dt)(di) 2 , which after an integration by parts 
gives — (l/2)[(d/dt)(g/g)]d 2 ; with g ~ e~^, the factor in square brackets becomes —0. Since 
= ^/Gl(X + Y), the rapidly varying X will average to zero, and one is left just with the 
very small Y, which is of order Hubble square. For the term (g/ g) 2 (a i ) 2 ) the term (X) 2 
also averages to order Hubble square, implying that the induced mass term is of horizon 
size. These "paraphotons" carry new relativistic degrees of freedom, which could in turn 
modify the Hubble expansion rate during Big Bang nucleosynthesis (BBN). Note, however, 
that these extremely light gauge bosons are thought to be created through inflaton decay 
and their interactions are only relevant at Planck-type energies. Since the quantum gravity 
era, all the paraphotons have been redshifting down without being subject to reheating, and 
consequently at BBN they only count for a fraction of an extra neutrino species in agreement 
with observations. 
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